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Abstract 

We consider ballistic aggregation equation for gases in which each particle is iden- 
tified either by its mass and impulsion or by its sole impulsion. For the constant 
aggregation rate we prove existence of self-similar solutions as well as convergence to 
the self-similarity for generic solutions. For some classes of mass and/or impulsion 
dependent rates we are also able to estimate the large time decay of some moments 
of generic solutions or to build some new classes of self-similar solutions. 



1 Introduction 

In the present work, we are concerned with ballistic aggregation Smoluchowski like models 
for which we establish quantitative information on the qualitative behavior of solutions. By 
ballistic aggregation, also (improperly) called kinetic coalescence in previous works [2j [7], 
we mean that we consider a system of particles identified by their mass and impulsion 
which undergo an aggregation mechanism. That differs from the simplest aggregation 
mechanism introduced by Smoluchowski [TS] in which model the particles are identified 
by their sole mass. 

Let us be more precise. We denote by P = P y with y = (m,p) a particle of mass 
m > and impulsion p £ if The space of particles states is then Y = M + x M. d and 
the velocity of the particle P y is v = p/m. We assume that at a microscopic level (the 
level of particles) the rate of collision of two particles P = P y and P' = P y i is a given 
nonnegative function a = a(y, y') and when these two particles collide they join to form 
one aggregated particle P" = P y " in such a way that the mechanism conserves total mass 
and total impulsion. In other words, the microscopic mechanism reads 

a>(y,y') 
y > ± y" > 
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with y" = (m",p") given by 



m = m + m , p = p + p. 



It is worth mentioning that the above reaction dissipates kinetic energy since that, denoting 
by = m" |v' : '| 2 /2 the kinetic energy of particle P", we have 

£**_£_£* = 1 \P+P*\ 2 1 \P\ 2 jWf 

2 m + m* 2 m 2 m* 

1 mm* 2 

2 m + m* 

At the mesoscopic (or statistical or mean field) level, the system is described at 
time t > by the density function f(t,y) > of particles with state y £ Y. For a 
given initial distribution the evolution of the density / is described by the Smolu- 
choswki/Boltzmann like equation: 

d t f = Q(f) in (0,+oo)xF, (1.1) 
/(0) = f m in Y. (1.2) 

The collision operator Q(f) is given by Q(f) = Qi(f) — Q2(f), where 

Qi(f)(y) = \ I r <y'i y - y') ftf) f(y - y') dm ' d p> ( L3 ) 



Q2(f)(v) = / / a(y,y') f(y) f(y')dm'dp'. (1.4) 
JR d JO 

The two following examples of functions a have been considered in relation with models 
in astrophysics [T71 18] : 

a(y,y , ) = a H s(y,y') ■= (m 1/3 + m' 1 / 3 ) 2 \v - v'\, (1.5) 

a(y, y ) = OArp(y, y ) := — , ^ . (1.6) 

mm \v — v \ z 

This model is seen as a simple test case or elementary analog of more realistic situations 
in fluid mechanics or astrophysics [H [9]. We refer to the introduction of |14} El [7] for an 
elementary introduction to physics motivation of such a model. We also refer to [U [9j [TBI 
[17] and to the references quoted in [TU [2J [7j for a more detailed discussion about physics 
of aggregation. 

In the context described above it is very natural to impose on the initial data fi n to 
have finite number of particles and momentum. This condition reads: 

< /* n € L 1 (y, (1 + m + dydp) . (1.7) 

Existence of solutions under that condition has been proved in [TU [21 [7] . It has also been 
proved that 

f(t,-)->0 in L l (Y), ast^+oo, (1.8) 
that is that the total number of particles tends to 0. 
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A more detailed description of the asymptotic behaviour of the solutions may be ob- 
tained by considering scaling invariance properties of the equations. This may be done 
for example by studying the so-called self similar solutions as it is possible to do for the 
Smoluchowski equation, see [H |6] and the references therein for recent results in that di- 
rection for that model. A first difficulty to this end is to determine the relevant scalings of 
the equation (jl.ip . (jl.3p . (jl.4p . We are very far from being able to treat the general case 
when the aggregation kernel a(y, y') actually depends on both mass and momentum of the 
two colliding particles and even when the aggregation kernel a(y, y') only depends on the 
momentum of the two colliding particles. We then may be less ambitious and just ask for 
whether a more accurate version than (|1.8p for some rate of aggregation a is available? 
We may imagine to answer that question in several ways listed below by order of accuracy, 
and indeed depending of the case we will establish any of such a kind of information. 

• Answer 1. Upper bound on moment: 3 a, 3 v, C 6 (0,oo) such that 

M a (/(t,.))<^ Vt>l. 

• Answer 2. Upper and lower bound on moments: 3 a, 3i/j = f«(a), Cj = Cj(a) (E 
(0, oo), such that 

^<M a (/(t,.))<§ Vt>l. 

• Answer 3. Existence of self-similar solution: there exists some profile function 
Voo : Y — * R+, some exponents A, fj,, v £ R such such that the function 

<p(t,m,p) :=t x <p 0O (t> J, m,t u p) 

is a solution to equation (jl.ip . (|l,3p . (jl.4p . 

• Answer 4. Convergence to self-similarity: for any given solution / there exists a 
self-similar solution ip such that / ~ (p as t — > oo, in a sense to be specified. 

Here depending of the model, we define the moment of order a of / in the following way: 

• when / = f(y) with y = m £ Y = (0, oo) or y = p G Y = M. d , then a = a£l and 

M^f) = M a (f) = j Y \y\ a f dy, (1.9) 

• when / = f(y) with y = (m,p) 6 7= (0, oo) x M d , then a = (a,f3) G R 2 and 

M & (f) = M a ,/3(f) = J Y m a \pf fdy. (1.10) 

The results obtained in this work are very partials and may be classified as follows. 

In Section 2 we consider the case of the kernel ausiUi y') (which depends on both mass 
and momentum) and the only result we are able to prove is a upper estimate on some 
moments (that is a result of type "Answer 1"). 
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In the remainder of the paper, we focus our attention on some toy models in which 
the aggregation rate a depends upon the only impulsion or upon the only masses, namely 
a {y^y') = a {PiP')i a {y^y') = a{ra^m!) or even a(y,y') = 1. The relation with the initial 
problem is not clear, and in particular it seems that a velocity depending aggregation 
rate a(y,y') = a(v,v') should be more natural that an impulsion depending aggregation 
rate a(y,y') = a(p,p'). Anyway, on the one hand such kind of aggregation rates has 
been considered by physicists, see [TJ [UJ [TBI E]) an d on the other hand our results and 
methods can give some ideas in order to tackle the so much more difficult models where 
the aggregation rate depends on both mass and momentum. 

Then, in Section 3 we consider a class of kernels which only depend on the momentum 
p and p' , we establish some moment estimates of type "Answer 2", from which we deduce 
the rather strange conclusion that solutions do not enjoin a self-similar property (nor self- 
similar solution exists). That result sow doubt about the fact that in the case of the mass 
and impulsion hard spheres kernel, solutions develop self-similar behavior. 

We treat in Section 4 the case where the kernel depends only on the masses m m! of 
the colliding particles and we exhibit a new class of self-similar solution (that is " Answer 
3"). Lastly, in Section 5 the case of constant kernel is treated, for which results of type 
"Answer 3" and "Answer 4" are established. 

We end that introduction by some remarks and open questions. A common feature of 
these equations is that 



and when the cross-section a is homogeneous of order 7 (which belongs to R or M?) it is 
likely that 



a result which is also known to be true for the coagulation equation (see El Hj) and for 
the inelastic Boltzmann equation (see and the references therein). The equivalence 
(II. lip is established for the the impulsion depending and the mass depending aggregation 
rate, but only one side of that equivalence is proved in the case of the true mass and 
impulsion depending hard spheres aggregation rate. We ask then. 

Open question 1. Is it true that the assymptotic equivalence behavior (jl.lip holds for 
some true mass and impulsion depending aggregation rate? 

An other interesting question should be to establish some asymptotic behavior of typical 
velocity or impulsion depending quantity. A way to express that in mathematical terms 
is the following: 

Open question 2. Is it possible to exhibit some moment for which we may determi- 
nate the long time behavior of M^/Mq (even just saying that it converges to something)? 

2 Mass & impulsion dependence remark on the hard 

spheres model. 

Let us recall the following result 



Mi, (t) = Mi, o (0) 



and Mo,o(i) — ► as t — » 00, 



My(t) = - as t — > 00, 



(1.11) 
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Theorem 2.1 [7| Theorem 2.6, Theorem 2.8 and Lemma 3.3] Assume that a satisfies 

< a{y,y) = a(y \y) < k s (y) k s (y'), Vy,y' G Y, 

a(m, —p, m , —p) = a(m,p, m ,p) V (m,p), (m ,p) G Y, 

a(m,p,m' ,p') < a(m,p,m , —p) V (m,p), (ra ,p) £ Y s.t. (p,p) > 0, 

with k${y) := 1 + m + \p\ + For any even (in the p variable) initial condition 

< fin G I jl (Y;kg(y)dy), there exists a unique solution f £ C([0, T); L (Y; k${y) dy)) n 
L°°(0, T; fc|(y) dy)) VT > ; which furthermore satisfies 



f{t,.)mdy= Cst, (2.1) 

f(t,.)is even, so that / f(t,.)pdy= 0, (2-2) 

/(t,.)l«| fc ^< [ fin\v\ k dy, Vfc>0, (2.3) 
y Jy 



f(t,.)\p\ 2 dy< / / in |p| 2 dy, (2.4) 
y iy 

/ fm a dy^0 when t — > oo, Va< 1. (2.5) 
iy 

Remark 2.2 (ij It is worth mentioning that the hard spheres collision rate ans does 
satisfy the assumption of Theorem \2.1\ but not the Manev rate onp- 
(ii) As a consequence of \2. \2. 3\) . \2.J$ and \2.5\) we deduce that 

M ai(3 (t) := J f(t, .)m a \pfdy as t oo (2.6) 

whenever (a, (3) belongs to the region 

{/3 G [0, 2], a < 1 - 13/2} U {(3 > 2, a < 2 - (3}. 

In the case of the hard spheres model we are able to quantify the rate of decay of one 
of the moment functions of the solution. More precisely, we have the following result. 

Lemma 2.3 Assume that a = aus- With the assumption of Theorem \2.1\ there holds 
A 1 := M_!/ 3)1 (0) < oo and 

M_ 1/3 ,i(t) < -x^m v ^°- ( 2 - 7 ) 

Proof of Lemma 12. 31 First we have M_ 1 / 31 (0) < oo because 

m _1/3 \p\ = m 2/3 \v\ < m 4/3 + \v\ 2 < 2k 2 s . 
Now, from the expression (|l.ip - (|1.2p of the collision kernel we have 

J y Q(f, f) m- x ^ \p\*V = \\ \ A_ 1/3> x f f dydy', 
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with 

A_ 1 / 3)1 = [(m + m')~ 1/3 Ip + p'l - m~ 1/3 |p| - (m')~ 1/3 [r + r'] 2 \v - v'\. 
On one hand — A_i/ 3il > because 

(m + m ' )V3 (JJs + (Jyw) * N + * |p + p/| • 

On the other hand, if we only take into account the values of v and v 1 where v-v 1 < and 
suppose that, for example, \p\ = min(|p|, \p'\) we have 

- A -A' s (^ + (^-(^T^))l'' 2 + <'''» 2 llH + Ml 
> ( M.) tt/f] tun = M M_. 

ym 1 / 3 / m 1 / 3 (m') 1 / 3 

Whence, by evenness of / 

d f f \p\ a s 1 f \p\ f t>A a i 

<- -K//7^*) 2 ' 

from which (|2.7p straightforwardly follows. □ 

3 The impulsion dependence case a = a(p,p*) 

We consider now the equation (jl.ip . (|1.3|) . ()1.4p with a collision kernel a independent of 
the mass of the colliding particles. We may then integrate the equation with respect to the 
mass and obtain that the function of t and p, J* °° fit, m,p) dm, that we shall still denote 
/, satisfies the equation: 

dtf = Q(f,f) m (0,+oo)xM d , (3.1) 
/(0) = f in in R d , (3.2) 

the collision operator Q{f) is given by Q(f, f) = Qi(f, f) - Q 2 (/, /), where 

QlifJM = 11 a(p',p-p') f(p') f{p-p')dp', (3.3) 



2 . 

Q 2 (/,/)(p) = / a(p,p') f{p) f(p')dp'. (3.4) 
We focus on the cases 

a(p,p') = \p-p'p, 7 €[0,2], d £ N* . (3.5) 

Before stating our main result we need some definitions and notations. We say that a 
function / on M. d is even if 

f(-p) = fip) Vp G R d , 
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it is radially symmetric if 



f(Rp) = f(p) Vp G R d , R G SO(d) 

where SO(d) stands for the rotation group on M. d . For any weight function k : M. d —>■ R+ 
we define the "moment of order k" of the non negative density measure/ G M} (R ) by 

M k (f):= [ k(p)f(dp), 
JR d 

and we define M\ as the set of Radon measures fi such that Mk(\fi\) < 00. For any a G R + 
we use the shorthand notation 



M a := / f(p)\p\ a dp, 



that is M a = M k (f) for k(p) = \p\ a and the shorthand notation M\ = M\ for i(p) 
1 + \p\ a . 



Theorem 3.1 Consider the aggregation rate f5T 

(i) For any even initial datum fi n G M\ a , a £ N\{0, 1}, there exists a unique even 
solution f e C([0,T);M 1 (R d ) - weak) n L°°(0, T; M 2 1 a (R d )) to equation if^7])-(CT. For 
any a G [0, 1] i/ie function t 1— > M a (t) is decreasing and f(t, .) is radially symmetric for 
any t > if furthermore fi n is radially symmetric. 

(ii) Moreover, the solution f{t,.) satisfies 

< M y (t) < ^ * ■ - - Vt>0, (3.6) 



M 7 (0)~ 1 + kit ~ 7W ~ M 7 (0)" 1 +M 
/or some constants ki = ki{^,d) G (0, 00). 

One of the main tools in order to establish that result is to consider moment equations. 
As it is classical for the coagulation equation, but here using one more change of variable 
p' — > —p', any even solution / to equation (|3.ip - (|3.4p satisfies (at least formally) the 
fundamental moment equation 



d 1 

dt. 

1 
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M a = \ I [ ff'a(p,p')[\p + p'\ a -\p\ a -\p'\ a ]dpdp' 

ff Hp,p') [\p + pT-\p\ a -\pT] 

+ a(p, -p') [\p - p'\ a ~ \P\° ~ \PT) } dpdp'. (3.7) 



More precisely, we consider in this Section the case 7 G (0,2) and d £ N*, the case 
7 = 1 and d = 1 and the case 7 = 2 and d G N*. The case 7 = and d = 1 is treated in 
Section 5. We shall use the following notation for the moments of order a G N: 
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3.1 Proof of the existence and uniqueness part in Theorem 13.11 

We prove in this subsection an uniqueness and existence result for a general class of 
aggregation rates by adapting some arguments from [10j Ej , see also [13] . We then deduce 
the existence and uniqueness part in Theorem 13. 1[ 

Lemma 3.2 We consider a continuous aggregation rate a : R 2d — ► R+ which satisfies 

a(-p,-p') = a(p,p') Vp,p'ER d , (3.8) 
a(p,p') < a(-p,p') Vp,p' G R d , p ■ p' > 0, (3.9) 

a even weight function k : R d — > R + and we define 

&k(p,p) ■= a(p,p) [k(p") + k(p) - k(p)}, A k (p,p') = A(p,p) + A(-p,p). 

We assume that 

a(p,p) <Ck(p)k(p) and A(p,p')<Ck(p)k(p') 2 . (3.10) 

For any given even initial datum fi n G M^(M. d ) there exists no more than one even solution 
f G C{[0,T);Ml(M. d ))nL oo (0,T;Ml 2 {R d )) to equations 



Remark 3.3 (i) The same result holds without the evenness assumption on the density 
function when the second condition in \3.10\) is replaced by 



/\2 



A(p,p ) < C k(p) k(p' 

We refer to \1(K [?]/ where such kind of result is proved in a L 1 framework. The same result 
also holds for radially symmetric solutions when we assume that 

a(Rp,Rp') = a{p,p) \/p,p' G R d , R G SO(d), (3.11) 

and the second condition in \3. 1 0\) is replaced by 

A(p,Rp')dR < Ck{p)k{pf. 



L 



ReSO(d) 

(ii) The same kind of result holds for aggregation rate defined on Y 2 with Y = (0, oo) x IR^ 
as it is the case when particles are identified by their mass and impulsion, see [7|/. 

Proof of Lemma [331 Step 1. We claim that for g G C([0,T);M^ - weak), G G 
L^OjTiMjJ) and b G C((0,T) x R d ;R + ) such that 

d t g = G-bg in the sense of V([0, T) x R d ), (3.12) 

the differential inequality 

^\\gk\\Mi < \\Gk\lMi - \\bgk\\ M i (3.13) 

holds in the sense of T>'([0,T)). First, it is clear using a classical duality argument 
that equation (|3.12p has at most one solution. Indeed, given two solutions 51,52 G 
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C([0,T);Mj£ — weak), we have for any t G (0, T), G C comp (R d ) and denoting by 
if G C comp ([0, t] x IR rf ) the solution to the dual homogeneous equation dt<p = bip 



(92 ~ gi)(t) (ft dp = / / {(d s g 2 - d s gi)ip+ (g 2 - gi)d s ip}dsdp = 0. 

JO JR d 

Now, for any g £ (0) G Cjc e ■= {u G C(R d ); suppn C K £ }, with K £ C M. d a compact, and 
any G £ G L 1 (0, T; C^) there exists a (unique) solution g e G C([0, T); Cr-J to equation 
(|3,12p which furthermore satisfies 



— [ \9e\kdy = [ (G £ - b g £ ) signg £ k dy 

JR d 



dt 



< [ \G £ \kdy- [ \g £ \akdy. (3.14) 

jR d JR d 



Here, signg £ = 1 if g e > 0, signg £ = if g e = 0, signg £ = — 1 if <? £ < 0. Finally, we can build 
(by a standard truncation and regularization by convolution process) the sequences (G £ ) 
and g £ (0) such that furthermore G £ — 1 G, g £ (0) — 1 g(0) in the weak sense of measures in 
Ml ||G e ( S )|| A/ i < \\G(s)\\ M i for a.e. s G (0,T), || 5e (0)|| M i < ||ff(0)|| M i. By the previous 
uniqueness argument we have g £ — g in the weak sense of measure and we get (j3. 13|) by 
passing to the limit in (|3.14p . 

Step 2. Let us consider two solutions fi,f 2 G C([0, T); M^(R d )) n L°°(0, T; M^ 2 (lR rf )) 
which are evens and let us denote -D = f 2 — fi, S = fi + f 2 - By a standard algebraic 
computation D satisfies the following equation 

d t D = Q(f 2 ,f 2 )-Q(fi,fi) = Q(D,S) 
= Q 1 (D,S)-SL{D)-L(S)D, 

where 

Qi((p,i>) = -(Qi(tp,ip) +Qi(ip,<p)), L(<p) := a(p,p')<p(p')dp'. 

1 JR d 

Because of the assumption made on a and / we have D G C([0,T];M^ - weak), G := 
Qi(D,S) - SL(D) G L°°(0,T;Ml) and < b := L(S) G C([0,T] x R d ) so that the first 
step implies 

i\\D\\ Ml < \\(Q 1 (D,S)-SL(D))k\\ M1 -\\DkL(S)\\ M1 



dt 



\jj a [k" + k']\D{dp)\S{dp')-\ j j ak\D(dp)\S(dp') 



< 

~ 2 

< \JJ A\D(dp)\S(dp') <j\\S\\ M i 2 \\D\\ M i. 

Uniqueness follows by using the Gronwall lemma. □ 

Lemma 3.4 Consider a continuous aggregation rate a : M. 2d — ► K+ which satisfies \3. 8\) 
(resp. \3.11\) ). i fPj) as well as 

a{p,p')<C(k + k') Vp, p G R d , (3.15) 
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for the weight function k(p) = 1 + \p\ 2 and some constant C G (0, oo). For any given 
even (resp. radially symmetric) initial datum fi n G M2 a (M rf ) there exists at least one even 
(resp. radially symmetric) solution f G C([0, T);M 1 (R d ) — weak) n L°°(0, T; M2 a (R d )) to 
equation \3. l\) - (3^\ ), and this one furthermore satisfies t i— ► Mg(i) zs decreasing for any 
/?€ [0,1]. 

Remark 3.5 is ZifceZy i/iai 6y adapting some arguments introduced in fl*$ . see also |3j, 
/or any even (resp. radially symmetric) initial datum /j n G L2a(^ rf ) ^ e approximating 
solution f n (t, .) 6mZi in i/ie proof below is a Cauchy sequence in C([0, T; L 1 (M d )) so i/iaf 
we may conclude f G C([0, T); L 1 ^)) n L°°(0, T; L^ a ( 



Proof of Lemma 13.41 We define the sequence of bounded aggregation rates a n := a An, 
for which classically fixed point argument (see for instance [7] which deals with some 
similar situation) implies the existence of a unique even (resp. radially symmetric) solution 
f n G C([0, T); L2 Q (lR d )) to equation (|3.ip - (|3.4p associated with a n for any initial datum 
fin,n G Ll a+2 (R d ), a G N, a > 2. Then, we have for any /3 G N*, j3 < a 



|//n(l + N 2/? ) = \ffnf n an [(\p\ 2 + 2p-p'+\p'\ 2 f-\p\^-\p'\^- 
= j fnf>n [2 (3 p-p'\ P \W-V -1/2 

where in the last sum we have Pi + P2 + P3 = and (Pi > 2 or (P2 > 1 and > 1)) or, 
in other words, \p-p'\^ |p| 2/?2 |p| 2/3a < |p| 2/3 \p'\ 2 (P~@ ) with 1 < P' < P — 1. Since we also 
have 

fnfn( a (P,p') An - a(-p,p') A n) p ■ p \p\ 2{f3 ~ l) < 0, 



p-p'>0 



we conclude with 
~dt 



J \<&<B-\ J 



a\p\ 2 ? \p'\ 2 ^'\ (3.16) 



l</3'</3- 

When P = 1 the set of admissible values of /?' is empty, and we recover a result from [lj 



tft / /n( 



l + |p| 2 ) < 0, 



so that 

sup \\fn\\ L \ < ll/m,nlLi- (3.17) 
[0,T] 

When P >2, gathering (|3. 15|) . (|3.16)) and ()3.17p . we easily conclude by a iterative argument 
that 

SU P ||/ n |U <C T (P,\\fin,n\\L\ )• (3-18) 
[0,T] ^ fc/3 
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Considering a sequence (fi n ,n) such that Un,n fin in the weak sense of measure and 
H/innllL 1 remains bounded, we easily pass to the limit in the equation satisfied by f n 

thanks to (|3.18p . The fact that t \—> Ma{t) is decreasing comes from the fact that p i— > \p\@ 
is a sub-additive function when (3 G [0, 1], so that A^ < and then d/dtMp(t) < 0. □ 

Proof of the existence and uniqueness part in Theorem 13.11 It is clear that 
a (p,p') = \p~ p'V satisfies (j3.8|) . (|3.9p . the first inequality in (|3.10p and (j3. 15j) . Moreover, 
the second inequality in (|3.1U|) holds since we have 

A 2 = \p-pA 1 (\p + p*\ 2 + \p*\ 2 - \p\ 2 + 1) + \p + p*V {\p-p*\ 2 + \p*\ 2 - \p\ 2 + 1) 

= 2 (\p -p*| 7 - \p + p*\' y )p -p + (|p - £>*| 7 + |p + p*| 7 ) (2 |p*| 2 + 1), 

where the first term in non positive and the second term is bounded by say 8 (k') 2 k, using 
that \p ± p*| 7 < 2 (|p| 7 + |p'| 7 ). We conclude by using Lemma 1331 and Lemma f3~4l □ 

3.2 Proof of the rate decay part in Theorem 13.11 when 7 < 2. 

For an even initial datum /j„ G Ml(WL d ) we consider the unique even solution / G 
C([0,T);M 1 - weak) n L°°(0, T; M\), VT, given by Theorem O^i) . This one satisfies 
the moment equation 

-^M 7 = i/ / ff^dpdp' = - f [ ff'A^dpdp', (3.19) 

with 

-A 7 = |p — p'| 7 [|p + j/| 7 - |p| 7 - |p'| 7 ] 

and 

-A 7 = |p-p'| 7 [Ip + p'I 7 - H 7 - b'l 7 ] +\p+p'V [\p-pV- W- b'PI • (3-20) 

We split the proof of Theorem I3.1( ii) in several steps. 

Step 1. One the one hand, for any given A > and any p,p' G M rf such that A' 1 \p'\ < 
\p\ < ^4 |p| we easily get 

|A 7 | < (\p\ + \p'\y max [(\p\ + \p'\y, |p| 7 + |p'| 7 ] 

< 2 4 max(|p|,|p'|) 27 < 2 4 A 7 (|p| |p'|) 7 . (3.21) 

On the other hand, we define M := max(|p|, \p'\), m := min(|p|, x := m/M G [0, 1], 
e := p ■ p' £ [—1, 1] and we compute (in the first line we have assumed that \p\ = M which 
is not a restriction to the generality because of the symmetry of A~) 

— A 7 = M 27 { \p — xp'Y 1 [l + x 7 — \p + xj5'| 7 ] + \p + xj5'| 7 [l + x 7 — \p — xp'| 7 ] } 
= M 27 {(l + x 7 )[(l + 2ex + x 2 ) 7/2 + (l-2ex + x 2 ) 7/2 ] 

-2 (1 + 2ex + x 2 ) 7/2 (1 - lex + x 2 ) 7/2 } 
= M 27 (2x 7 + 0(x 2 )} < 3M 27 x 7 = 3(|p| |p'|) 7 (3.22) 

uniformly on e G [—1, 1] and x < A^ 1 for > 1 large enough. 
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Gathering (pT2TJ) and (pT22j) we obtain 

±A 7 >-h \p\-r\pV 



\/p,p' G 



with k\ := max(3/4, 2 3 Aq)/4, and equation (|3.19p then implies 

jM„, > -k x M 2 . 

We straightforwardly obtain the first inequality in (|3.6p by integrating this differential 
equation. 

Step 2. First, together with the variables M, x and e introduced in Step 1, we define r > 
and u G [0, 1] by setting r 2 := \p\ 2 + |p'| 2 and u := 2p -p' /r 2 , so that |p±p'| 2 = r 2 (1 ± u). 
Splitting the positive and the negative terms in identity (|3.20p . we have 

-A 7 = W + \p'[ 1 ){\p-p'\ 1 + \p + p'[ 1 )-2\p-p'[<\p + p'[< 

( (\p\ 2 p /2 + (\ P fr /2 



,27 



(l+«)^ + (l 



\ (H 2 + |p'l 2 ) 7/2 

Since 7/2 G [0, 1], the map x 1— » x 7- / 2 is sub-additive, and we obtain 



2(l+u) 7 / 2 (l 



-A- 



> r 



2t 



{ [(1 + uy/ 2 + (1 - u) 7 / 2 ] - 2 (1 + u) 7 / 2 (1 - u) 7 / 2 } 



> M 27 (l + -u) 7/2 (l-n) 7 / 2 0(-u), 0(«):= (1 - u)^'' 2 + (1 + u)^'/ 2 



We easily verify that (j) is increasing on [0, 1] so that 4>{u) > 0(0) = for any u G [—1,1], 
Coming back to the variables M, x and e, that is (j)(u) > for any p,p' G M rf such 
that the associated variables M, x and e satisfy M > 0, x > and e / 0. Moreover, when 
e = (p and are orthogonal vectors) we also have 



-A, 



+ \ P fr /2 



\pV + \pT 



+ Ip?)^ 



> 2M 27 



1 + x 7 - (1 + x 



2n 7 /2 



> 



for any p,p' G M. d such that the associated variables M and x satisfy M > 0, x > 0, because 
the function z z" 1 ' 2 is strictly sub-additive on R+, that is (z + z')^l 2 < z~<l 2 + (z') 7/2 
for any z, z' > 0. Gathering these two lower bounds on — A 7 , it yields 

- A 7 > M 27 V(x,e) (3.23) 

with i/j(x, e) > for any x > and e G [—1, 1]. 
Next, coming back to (|3.22|) . we also deduce 

- A 7 = M 2j {2x 7 + 0(x 2 )} > M 27 x 7 (3.24) 

uniformly on £ 6 [-1,1] and x < A^ 1 for Aq > 1 large enough. Gathering (|3.23p with 
(|3.24p we deduce that for some constant k% > we have 



Vp,p G R d 
and equation (|3.19p then implies 



^A 7 > k 2 M 2 ^x^ = k 2 (\p\ \p'\r, 



d 



Mv < -koM 2 



dt 7 ^ 

The second inequality in (|3.6p is again obtained by integrating this differential equation. 
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3.3 The case a(y,y') = \p — p'\, d = 1. 

In the particular case under consideration d = 1 and 7 = 1, we can establish a more 
accurate version of the decay estimate on the moment M\ together with additional moment 
estimates. 

Lemma 3.6 Assume a(y,y') = \p — p'\ and d = 1. For any even initial datum fi n G 
M3 (R) the unique solution f G C([0, T]; M 1 (M)) n L°°(0, T; M3 (M)) 0/ (E2P-|53P #roen 
fry Theorem \3.1\ satisfies for any t > 

f Mq(0) 2 3 / 2 M (0) \ M (0) 

maX l v (l + M 1 (0)W (2 + 3M i/3 (0)t) 3/ 2 J * MoW * (l + M l( 0)t)V 2 ^ 

— — i < Mi(t) < — \ (3.26) 

A/2(()) < M 2 (t) < M 2 (0) (3.27) 

< M 3 (t) < M 3 (0). (3.28) 



(l + Mi(0)t/2) 2 

M 3 (0) 
(l + M!(0)t/2) 2 



Remark 3.7 77ie above estimates on the behaviour of M\(t) for t large are quite good. 
That is not the case for the estimates on M a , a = 0,2,3 which seem to be rather partial. 
Worst, with these bounds we can not even know what is the limit of any of the quotients of 
moments M a (t)/Mi(t) for a = 0, 2, 3 as t — > 00. The value of such a limit would indicate 
whether the solution f(t) has a tendency to concentrate or to spread as t increases (see 
also below the discussion concerning the case 7 = 2). 

Proof of Lemma 13.61 Introducing the notations M = max(|p|, \p'\), m = min(|p|, \p'\), 
we systematically exploit the differential equation 



Jt M " = 4 



\ [ I ff'A a dpdp> (3.29) 
4 JrJr 



with 

A a ■= [M - m] \{M + m) a - M a - m a ] + [M + m] [(M - m) a - M a - m 



Step 1. a = 1. We have 

Ai = -2(M + m)m, 

from which we deduce 

|M 1 (t) = -^M-^, Bl (t):= ! [ ff'{min(\p\,\p>\)} 2 d P d P >. 

Since < {min(|p|, |p'|)} 2 < \p\ \p'\, we have < B\{t) < Mf(t) and we obtain the two 
closed differential inequalities 

-M 2 (t)<|M l( t)< J^M, (3.30) 
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from which we deduce (|3.26p . 
Step 2. a = 0. We have 

A = -2M, 

from which we deduce 

±M (t) = -^P-, B Q (t):= [ [ ff max(\p\,\p'\)dpdp'. (3.31) 
at 1 JrJr 

Since \p\ < max(|p|, \p'\ ) < \p\ + \p'\, we have M M x < B < 2M Mi and then 

- M Mi < -^M < M Mi. (3.32) 

Using the previous estimate (|3.26p on M\(t) we get 

M (t) < d < M (t) 



Mf x (0) +t/2 ~ dt ~ 2(Mf 1 (0) +£)' 

from which we deduce the first lower estimate as well as the upper bound in (|3.25p . 
Step 3. a = 2. We have 

A 2 = -4mM 2 



from which we deduce 



j t M 2 (t) = -B 2 (t), B 2 (t) :-- 



/ / //' min(|p|, \p'\) \p\ \p'\dpdp'. 
JrJr 



Using that < min(|p|, \p'\) \p\ \p'\ < \p\ 2 \p'\ together with (|3.26p . we obtain 

-M 2 — — ^ t < -M% Mi < — M 2 (t) < 0, 

Mi(O)" 1 +t/2 ~ ~ dt Iy ' ~ ' 

which implies (|3.27p . 
Step 4. a = 3. We have 

> A 3 = -2Mm 3 -2m 4 > -4Mm 3 > -4\p\ 3 \p'\, 

from what we deduce 

0> j t M 3 (t) > —Mi M 3 , 

which again implies (|3.28p . 

Step 5. a = again. Coming back to the moment Mq, we write for any e > 



= ~\ i i f f'lP' -P\dpdp' 
z JrJr 

- ~7 I I f {£ + -\p-p'\ 2 )dpdp 
4 JrJr £ 

> - S -Ml--^M Q M 2 . 
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By interpolation we have M 2 (t) < M Q 1/3 (t) M 2/3 (t). Since, by M 3 (t) < M 3 (0) for 

all t > we deduce M 2 (t) < M Q 1/3 (t) M 2/3 (0). Therefore 

jM (t) > - £ -M>-Im^M^(0) 

We now chose e = e(t) > such that eM 2 = ^ M^ 3 M^ 3 (0), or equivalently e = 

1/3 i /3 

Mn M 3 (0). With that choice of e(t) the moment equation reads 

| Mo (t)>-^M 3 1/3 (0)M 5/3 , 
from which we deduce the second lower estimate in (|3.25|) . □ 

Remark 3.8 In the last step, we may also argue as follows. Gathering the estimate 
m&x(\p\,\p'\) > (\p\ Ip'l) 1 / 2 , the differential equation h3.31\) and the interpolation estimate 
Ml 12 < Mf, 2 All 12 we obtain thanks to (E 



Together with h3. 26\) we recover the second lower estimate in {2j 



3.4 The case a = \p — p* | 2 

In the particular case under consideration 7 = 2 and rf 6 N*, we can close the family of 
moment equations for any moments M 2a , ce E N. In the following lemma we give the 
expression of moments up to order 4, showing a (unexpected?) non self-similar behavior 
of solutions. 

Lemma 3.9 Assume a(y,y') = \p — p'\ 2 and d E N*. There exists a numerical constant 
kd E (0,oo), k\ := 2, such that for any radially symmetric initial datum /j n E Mg(M) the 
unique radially symmetric solution f E C([0, T]; M 1 (R)) n L°°(0, T; Mg (R)) of nn\)-(3Zfi) 
given by Theorem \3.1\ satisfies for any t > 

Ms «> = MitorW (3 - 34) 

M A (t) = M 4 (0)(M 2 (0)' 1 + 2k d t) 1/kd ~ 2 . (3.35) 
Proof of Lemma 13.91 We proceed in several steps. 

Step 1. a = 2. Using the fact that / is radially symmetric (so that the odd moments of 
/ vanish) and the notations p = ro~, r = \p\, p' = r' a' , r 1 = \p'\, the fundamental moment 
identity (|3.T() implies 

T + M * = \ I I ff'[\p\ 2 -1p-p' + \p'\ 2 \{1p-p')dpdpi 

at z JR d JR d 

f f [P-P? dpdp' 

f'OQ f'OO p p 

-2/ / f{r)f{r')r d+l {r') d+1 drdr'x / [a ■ a 1 } 2 dado' 

Jo Jo Js^Js 11 - 1 

-2k d Ml 




15 



with 



/ / [a- a'} 2 dada') x meas(S a 
Js d - 1 Js d - 1 J 

meas^- 1 )- 1 / of da. 



i S d-i 

We compute k\ = 1, k 2 = 1/2. The expression (|3.34p immediately follows by integrating 
that ODE. 

Step 2. a = 0. When a = 0, the fundamental moment identity (|3.T|) and the fact that / 
is radially symmetric imply 

|m = ~ // ff'[\p\ 2 -2 P -p' + \p'\ 2 }(-l)dpdp' 

at z JR d JR d 

= -M 2 M . 

Integrating that ODE with the help of (ETMj) we get (ET33J) . 

Step 3. a = 4. When a = 4, the fundamental moment identity (|3.7p and the fact that / 
is radially symmetric imply 



\l I ff [\P\ 2 ~2p-p' + \p'\ 2 } [4 (p ■ p'f + 8 \p\ 2 (p ■p') + 2 \p\ 2 \p'\ 2 } dpdp' 
ff {[2 \p\ 2 } [4 (p ■ p') 2 + 2 \p\ 2 \p'\ 2 } - 16 \p\ 2 (p ■ p') 2 } dpdp' 




dt 2 
1 
2 

= 2// //' {N 4 |p , | 2 -2|p| 2 (p-p / ) 2 } dpdp' 
= (2- Ak d )M 2 M A . 

Integrating that ODE with the help of (|3T34l) we get (^BTj) . □ 

Remark 3.10 (i) On the one hand, the moment M a (g(t, .)) 0/ a s elf- similar function g 
of the form g(t,p) = t^G(t u p) satisfies 

M a (g(t,.))=C a t^ d+a >. 

On the other hand, when d = 1 we have k± = 1 so that the solution f of equation 113. l\) - (3~^\) 
satisfies 

M (f(t,.))~C' t^ 2 , M 2 (f(t,.))~C' 2 t-\ M 4 {f{t,.))~C' A r l . 

Since the long time behavior of these functions are incompatibles, there does not exist any 
self-similar solution with self-similar profile G 6 Mg (R). 

(ii) When d = 1, to make the ideas simpler, the moment M 2a satisfies the edo 



d / 2 \ / 2 

T M 2a = £ " M 2/3 M 2 ( Q+ i_ / 3) "E (2/9 + 1 J M ^+ 2 M2 ^-^' 

3=1 V J 8=0 V 



/3 

In particular, we find 



^M 6 = 3M 2 M 6 -5Ml 
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When M2(0) = 1/2 (for the sake of simplification again), the solution is 



M 6 (t)= (m 6 (0)-2M 4 (0) 2 + ^) (1 + t) 3 / 2 Vi>0, 

with Mq(0) — 2M4(0) 2 > (Holder inequality). The solutions of equation H3.1\) - pT4\ ) have 
a rather strange behavior since that 

M ~K t- 1/2 , #~^- 1/2 > #~^ 3/2 > 

Mq Mq 

In some sense, the behavior is in part comparable with the solutions of the inelastic Boltz- 
mann equation which energy (here the Mi moment) dissipates and in part comparable with 
the solutions to Smoluchoski equation which high moments rapidly increase. It is worth 
mentioning that here the "mean second moment" (that is M~2/Mq) tends to in the large 
time assymptotic. The opposite feature occurs for some models dealt in section^ 



4 The mass dependence case a = a(m, m*) 

Consider now the problem (|l.ip - (|1.4p where the kernel a(y, y') only depends on the masses 
of the particles, namely 

a{y,y') = a(m,m'), (4.1) 
and introduce the associated Smoluchowski equation 

dF 1 f m 

-(£, m) = — / F(t,m — m')F(t,m') a(m — ml \rn') dm' 



dt K ' ' 2 



F(t,m)F(t,m') a(m,m') dm' . (4-2) 



For any function ip G L 1 (M 3 ) we define the Fourier transform T and the inverse Fourier 
transform J 7 ^ 1 by 



^(rf) = (JFV)(»7) = / ^(P) e~ ip ' V dp, i>){p) = (2<T d / ^(rj) e 1 ^ dp. 

Theorem 4.1 For any smooth function a on M 3 homogeneous of degree 9 , 6 G (0,oo), 
and such that ip := F^ 1 (e~ a ^) > 0, and for any solution F = F(t,m) to the coagulation 
equation with coagulation kernel a(m,m'), the function f(t,m,p) defined by 

f(t,m,p) = m- 39 F(t,m)ip(^), (4.3) 



vr 



/,s- 



a solution of the equation lll.l]) . \1.3\) . |i.^[ ) for the same aggregation kernel. 



Remark 4.2 Theorem \4-l\ is not a general existence result of solutions to M.l\) . 111.3]) . 

Notice indeed that the initial data satisfied corresponding to these solutions are all 
of the form /(0,m,p) = m~ 3e Fi n (m) 92 (p/m 9 ^. An example of admissible function a is 
a(p) ■= \p\ 2 , so that that 6 = 1/2. 
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Proof of Theorem 14.11 We have to check that the function f(t,m,p) defined by (|4.3 
solves CUE]), (USD, (HU). We start with writing 

p \ dF 



U J -39 
— = m ip , 

at V rrr J ot 



F(t, m — m')F(t, mf) a(m — m , m) dm' 

F(t, m)F(t, m) a(m, m') dm' 



(4.4) 



On the one hand, using that 



^)^ = ^)(0) = e- a W = l, 

'R 3 

the last term in ()4.4p gives 

m~ 3S <p ( -^-k ) / F(t,m)F(t,m!) a(m,m') dm' = 
\m° J n 



= m- 3e tp(^) F(t,m) [ aim, m') Fit, m') I (m')~ 36 <p(-^) dp' 
\m y J J J R 3 m m 

f'OO r 

= f(t,m,p) / a(m,m) f(t,m',p')dp'. 
Jo JR 3 

On the other hand, let us define the function 

g(m,p) = m~ 39 ip(p/m ). 

Using the definition of 92 and the homogeneity of a, it satisfies for any < m' < m 

g(m,rj) = (p(m 8 rj) = exp(— a{m d ry)) = exp(— m a(rj)) 
= exp(— m a{rf)) exp(— {m — m ) a(rj)) 
= 9(m',7])g(m-m',r]), 

or coming back to the origin function 



(4.5) 



\m 6 



g(m,p) = / g(m ,p) g(m — m ,p — p) dp 1 . 
Jr 3 

Using that identity in the first (gain) term in (|4.4p . we get 

F(t, m — m')F(t, m) a(m — m , m) dm = 

g(m,p) / F(t,m — m')F(t,m') a(m — m ,m) dm' 
Jo 

F(t, m — m)gim — m ,p — p) Fit, m')g(m ,p) aim — m , m) dm dp 
f(t, m — ml ,p — p') fit, m' ,p') a(m — m! , m) dm' dp' . (4-6) 



We conclude that / satisfies [HJ , flL3]), (H31) by gathering (^4"|1 . and gS]). □ 

The previous Theorem is useful in order to prove the existence of self similar solutions 
for some kernels a(m, m') as it is seen in the following corollary. 
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Corollary 4.3 Suppose that a and 9 are as in Theorem \4-l\ Assume further that F is 
a self similar solution of the coagulation equation with coagulation kernel a(m,m'). Then 
the function f defined by M-Sty is a self similar solution of hl.l\) . hi. Sty . |73 



(j4tfp_ 



Proof of Corollary 14.31 The hypothesis on F means that for some functions v(t) 
and n(t) it may be written as: 

F(t,m) = z/(t)$(/i(t)m). 
Therefore / is a self-similar function since it may be written as 

f{t,m,p) = m- w v{t)MAt)rn)v(-^) 

= v{t)fi{tf e (fi(t)my se $(n(t)m)<p 

= v{t)^{tf e m(n{t)m^(t) e p) 

with *(M, P) = M~ 3d *(M) ip (P/M 6 ) . □ 

Remark 4.4 (i) Self similar solutions for equation hl.l\) . hi. Sty . \l-4\ ) had already been 
obtained in Q/. They correspond to the case 6 = 1/2 of the above Corollary. 

(ii) Self similar solutions of the coagulation equation are well known to exist for the cases 
a(m, m!) = 1, aim, m 1 ) = m + m' and a(m, m') = mm 1 . Their existence for several other 
kernels with homogenetity A < 1 have been proved in J2[/ and J2|/. In that last case, these 
self similar solutions are of the form: 

F(t, m) = — — . (4.7) 



We deduce under the assumption of the above Corollary that 

f(t,m, P )=t-rhm- 3e $>(^-) <p(Ai). (4.8) 



irr 



is a self similar solutions to equation hl.l\) . hi. Sty . Ji.^[ ) for the same kernel a(m,m'). A 
straightforward calculation yields 

P k (t) = / / \p\ k f(t,m,p)dmdp = t~~ / \P\ k <p(P)dP / M ke $(M) dM. (4.9) 
Jm d Jo Jw. d Jo 

As a consequence, we have Po — > 0, P\ — > and more generally Pu — » whenever k < 
but Pk/Po °° f or any k > and Pj. — > oo whenever k < The rough physics 

interpretation is that the total number of particle decreases, the total impulsion of the gas 
also decreases, but for instance the mean second moment P2/P0 tends to infinity in the 
large time asymptotic, which is the opposite behavior with respect to the one discussed in 
Remark \3.1(K Here, the behavior is quite similar with the bahavior of the solutions to 
Smoluchoski equation since the mean impulsion moment Pk/Po —> 00 f or an V k > 0. That 
makes again a difference with the model discussed in Remark \3.10\ 
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5 The constant case a = 1 



For the sake of simplicity, we restrict our study to the case d = 1. It is likely that it 
extends to higher dimension d G N*. 

Theorem 5.1 Suppose that the initial data fi n is even, regular and good decreasing prop- 
erties. Then il.2\) . il.3\) jl-4\ ) has a solution given by: 

f(t,m,p) = F- 1 (£- 1 F)(t,m,p) (5.1) 

F(t,(,0 = ry^ F 777^' ( 5 - 2 ) 

(H + (t/2)) 2 I 



FlfiXT) H + {t/2) 
with Hq := Mofi(fi n )~ 1 as defined in l\2.6\) . Furthermore, f satisfies 

2 H% rn An 2 

t 5 ' 2 f(t,tm,y/ip) ^<Poo(m,p :=^=4= = (5.3) 

in the weak sense of measure a(M 1 (Y), C C (Y)) as t — > +oo, where 

poo p 

A = H 2 (0) / / mf(0,m,p)dpdm, (5.4) 
Jo Jr 

h 2 (o) f°° r 

B = / / p 2 f(0,m,p)dpdm. (5.5) 



2 

Proof of Theorem 15.11 We first notice that the equation (jl.ip . (|1.3j) (|1.4p is now: 

lff m 

dtf(t,m,p) = - J J f(t,m — m,p — p)f(t,m,p)dmdp 

p poo 

-f(t,m,p) / / f(t,m',p)dm'dp'. (5-6) 

jR d JO 



2 Jm<i Jo 



This equation may be explicitly solved using Fourier transform with respect to p £ R and 
Laplace transform with respect to m > 0. Of course this needs the transform F of the 
function / to be defined. This has then to be checked once the expression of / is obtained. 
We thus define 

poo p 

F(t,(,0= / e~ mC e- ip tf(t,m,p)dpdm. (5.7) 
Jo JR d 



We then take formally Fourier and Laplace transforms in (|5.6p to obtain the Bernouilli 
equation: 

d t F(t, (, = \F 2 {t, C, - M (t) F(t, C, (5.8) 
M (t) = F(t,0,0). (5.9) 

We first notice, taking ( = £ = in (pTgj) . that M (t) satisfies ^M (t) = -\M$(t) from 
where 

M » ( " = B^m- (!U0) 
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Classical ODE integration methods lead that the solution of (|5.8p is the function F(t, £, £) 
given by (|5.2p . On the one hand, the function 1 1— > (Hq t/2) / (Hp+t/2) is strictly increasing 
with limit in infinity equal to Hq, so that for any 5 G (0, 1) there exists T G (0, oo) 



Vt G [0,T] 



Hp t/2 
H + t/2 



<flo(l-*), 



and on the other hand 



f'OQ P 

|F(0,C,OI< / / f(0,m,p)dmdp = H 1 . 

JO 



(5.11) 



(5.12) 



Gathering (|5.1ip and (|5. 12j) the fraction in the right hand side of (|5.2p is well defined for 
all t > 0. More precisely for any t G [0, T] 



1 



#o*/2 



F(0,££) H + (t/2) 
which implies 



> 



1 



Hpt/2 
H + (t/2) 



> tocoi -1 - #o(i - > * w-cor 1 



l^(*,C,OI < 



5(#o + V2) 



:i^(o,c,e)i- 



(5.13) 



As a consequence, any "good" decay and regularity properties of the initial data /j n ensure 
"good" decay and regularity properties of F(0, £, £). ^ is then possible to take the inverse 
Fourier and Laplace transforms of F(t,C,£) to define the function f(t,m,p). 
If one is interested in the behaviour of f(t,m,p) as t — ► oo it is a classical argument to 
consider the rescaled functions ip associated to / by the relation 



<p(t, M, P) := t 5/2 f(t, tM,ViP), 



so that 



Taking the Fourier and Laplace transform in both side yields 

F(t,c,0 = t-^tcVtO 

with 



tflg 



(Ho + (t/2)Y 



1 Hp t/2 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



Since we are interested in the long time behaviour of $(-,C>£) f° r an C an d £ fixed , we 
may write: 



1 



Hp t/2 
H + (t/2) 



1 



Hp + 



Hi 



Ho + [t/2) 
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and consider the auxiliary function 



^/ 2 )) 2 (^txy-^+(fi) 

' (5.18) 



We perform the following expansion up to the order o(l/t): 

#o = 7^— (0,0,0) + ^——(0,0,0) + -^-—^(0,0,0) +o - . (5.19) 



F(0,f,^) u t d( " ' ' ' y/i di " ' ' ' 2 t 3£ 2 v ' ' ' V* 
Since by hypothesis / is even with respect to p, we have 

— (0,0,0) = -i/ / f in (m,p)pdpdm = 



and then: 



We also have 



^fr< ' '°>=-F(aWf< ao ' >= - (5 ' 20) 

-— (0,0,0)=-/ p 2 f(0,m,p)dpdm, 
d£ Jo Jr 

which with the help of (|5.20p implies 



<l£? (0.0.0) = "f" (0,0,0)-^(0,0,0) + F-^(0,0,0)l — (0,0,0) 

f'OO r 

= H 2 (0) p 2 f(0,m,p)dpdm =2B. (5.21) 

Jo Jr 

Similarly, we compute 

OF f°° f 

— (0,0,0) = - / mf(0,m,p)dpdm, 

Cs jo Jm. 

which implies 

^ (0 ' ' 0) = -FTOf (0 '°' 0,= - 4 - <5 ' 22) 
Thanks to ([539]) . flSZDD , ([OT]) and (|5T22l . we deduce that (f5TT8|) reads now: 

4ff 2 (0) 

( *' C ' C) ~ {CA + eB + 2H 2 {0) + o(l)) 

from where 

Jim *(t, C ,0 = t hm *(t,C,0 = AC + Be y 2Hm =■ *oo(C,0. (5-23) 
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In order to come back to the original variables, we recall that from standard integral 
calculus for any C,V > 

_J_ r f e -mC e -i P ^- Cm ^ dpdm= 1 

(2tt)V2 y j R e d p dm c+ve+c 

from where choosing C := 2 Hq/A and T> := B/A, we obtain 

|„|2 

(^- 1 £- 1 )(^ 00 )= - ° - _— =i Poo (m,p) 



(2tt) 1 /2^ 



as defined in (|5.3p . Finally, (|5.24p implies that .) — ipoo in the weak sens of measure, 
which is nothing but (|5,3p . □ 

The previous Theorem shows the convergence of some of the solutions of (|l.ip - (|1.4p to 
a function which is a self similar solution of the equation (|l.ip . (|1.3|) . (|1.4p . i.e. a solution 
of the form 

f(t,m,p)=t- a (p(t- 1 m,t- l3 p) (5.24) 

for some function (p. The numbers a and (5 define the scaling of the self similar solutions. 
In the Theorem 15. II we have a = 5/2 and (5 = 1/2. It turns out that equation (jl.ip . (jl.3p . 
(|1.4|) has more than one self similar solution with the same scaling as it is shown in the 
next Theorem. 

Theorem 5.2 Let $ G C 1 ^) is smc/i f/wzf 



3(y,a;) = 1 £ ( 1 



2C$(^)+1 



is well defined for x £ R rf and y > 0. Then 



(5,5) 

is a se// similar solution to U.l\) , il.3\) . 

Proof of Theorem 15.21 We look after self similar solutions of the form (|5.25p . The 
function g must then solve: 

5 1 1 f f y 

- 2 9 ~ ydy9 ~ ijxd x g = - \ \ g(y-y',x- x') g(y',y')dy' dx - 



o 



oo 

-all g(y',x')dy' dx'. (5.26) 
Jr Jo 

We integrate this equation with respect to x and y and obtain 

> 

g(y',x')dy'dx' = 2. (5.27) 
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We now Fourier transform with respect to x and Laplace transform with respect to y: 

C<%? + = \g 2 ~ 9- 

We divide by g 2 and define G = l/g: 

Cd c G + £dtG = G-± 
The function G may then be any function of the form: 



c 

for any arbitrary derivable function Therefore 



<K(,0 = A\ — > ( 5 - 28 ) 

+ 1 



with, due to (|527|1 : 

2 /£ 2 

lim t — r = 2 <^=^ lim — 

C-o, 2 c $ + ! C-o, f->o V C 



If we want to define the function g from (|5.28p the function <3? must be such that g has an 
inverse Fourier and Laplace transform. □ 

Remark 5.3 If ®(z) = z + l, 

9 9 

9(C,0 = 




This is the profile of the self similar solution which appears in Theorem I5.il It is easy to 
obtain particular solutions g, some of them are explicit others are not. If, for example, 
<3? = 1 then g(y,x) = e~ y 5 x= q. Another explicit example is for &(z) = z which gives 

g(x,y) = \pnb y= Qe v3. On the other hand, if we take $(z) = y/z, the inverse Laplace 
transform, let us call it h(y,£), is still explicit: 

( )= c Wc? 7 J? ' 

It remains to check that h(y, •) has an inverse Fourier transform with respect to the variable 
£. It is easily checked that, for all y > fixed: 

%'0 = ^(-^). as^O 

^ +o(y 
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This function is then in L 2 (R) with respect to the £ variable and has then an inverse 
Fourier transform with respect to £ which is g(y, x): 

g(y,x) =J c "^ 1 (/i(y,-))(x). 

Moreover, for all y > 0, g(y, ■) G L 2 (IR) and the convolution of g(y,-) with itself is well 
defined 

t (g(y - y', ■) * sOA ■)) (0 = % - iA 

and 

f ^ - y', •) * g(y', •)) (£)*/ = T % - 2/', Kv , dj/. 

JO JO 

Therefore, 

/ |^( ff (j/-j/,-)*^,-))(0|de < /V |%-y',6%',0|^^ = ^4, 

JR JO JR k=1 



with 



rv/ 2 r 

, , \h(y-y',0Hy',0\dCdy', 

ry/2 



h:= I \h(y-y',0Hy',0\dCdy', 

Jo V /2 <l5l<Q/-s/') 1/2 
rv/ 2 r 

h:= I i/9 \h(y-y',OHy',0\d^dy', 

Jo Jy^ 2 <(y~y') 1/2 <\e\ 



h := „ t \h(y-y',Oh(y',0\dtdy', 

Jy/ 2 J\S\<(y-y') 1/2 <y' 1/2 

h:= f V [ i/9 , \h(y-y',Oh(y',0\d^dy', 

Jy/2J( y - y >) 1/2 <\t\<y'^ 2 

h-= f V [ , \h(y-y',t)h(y',0\dtidy'. 
Jy/2 J(y-y') 1/2 <y' 1/2 <\(\ 

We must verify that each term is finite. Indeed, we have 

rv/2 r & 

h < C / — ; — 5 ^didy' 

Jo J\e\<v' 1 ' 2 <(v-v') 1 / a y ,3/2 (y - v'Y' 2 



'\i\<y' 1/2 <{y-y') 1 ' 2 y' 3/2 (y 
y/ 2 min{y' 3 ' 2 , (y - y'? /2 } , , 



h < C 



151 b^-^W 



ry/2 



o Jy> 1/2 <m<(y-y') 1/2 (y-y') 3/2 

C [ y ' 2 f /ll 



- Wi S <»-»r'AW + lii + l ^ /:;; " 
c r /2 r ( 2 

- 



, , j v I , {^n + x ) d ^ d y < °° ; 

y 3/2 Jo Jy'^m^y-v 1 ) 1 ' 2 VVy 7 ' 
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h - ° C 1<"< { ,-,;"><KI (WW\ + ? + I.? 1 1 x 



[ [ ( 1 1 ( 1 1 \ 

? + °(w) + (^) + (^))« 



C [ y / 2 If df 
— — / ~~f= / i — To dy + 

\/yJo vV A ,1/2 <(j/-j/') 1/2 <i?i ISI 



+C< i) ( + Sy' 1/2 <{y-y') 1/2 <\(,\ l^l 3 ^ + 



Jo A' 1 



_c_ r y / 2 d y r y i 2 f j_ j_\ z^ 2 

" vWo y/^Viy-y 1 ) + h \Vv + Vv 7 ) y-y ,+ Jo (y-y') 2 ' 

Similar estimates show that the integrals and Iq converge. The function Jq h(y 

2/'>£) h(y' \0 dy is then in L 1 (IR) and has then an inverse Fourier transform which is 



f\g{y-y',-)*g{y', ■))(£) 

JO 
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